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Vibration Compensation for Sensors 

Field of the Invention 

The present invention relates to sensors, and in particular to the 
5 compensation for vibration-induced error affecting the operation of such 
sensors. 

Background of the Invention 

IjpaaLsensors, such as gyroscopes and accelerometers, are often 
1 0 adversely affected by V1 brations. They may be very sensitive to such vibrations. 
In fact, vibration sensitivity can often be a dominant performance limiting error 
source. Many inertial sensors rely on micro-electromechanical sensor (MEMS) 
devices which utilize vibratory mechanisms for rate/acceleration sensing. These 
types of sensors have shown evidence of significant vibration sensitivity 
15 problems. 

Inertial sensors are often used in environments that inherently subject the 
sensors to significant vibrations. These include use in aircraft and weapons, 
where the power plants and environment serve as significant sources of 
vibration. For example, a commonly used resonant beam accelerometer that is 
20 often viewed as a one milli-g accuracy sensor can exhibit bias shifts on the order 
of several milli-g' s under typical operational vibration environments. 

For inertial sensors which exhibit vibration sensitivity, the most common 
effect is a slowly varying low frequency error component which changes as a 
function of the applied vibration spectrum. In some applications, such as a 
25 global position sensor aided inertial system, the effect of the vibration induced 
bias shift can be trimmed or estimated such as by use of a Kalman filter. 
However, when the vibration level or spectrum shape changes significantly, the 
vibration-induced error also changes, diminishing the value of any previously 
obtained calibration. 

30 In the case of an inertially guided munition launched from an aircraft 

equipped with a high quality inertial navigation system (INS), an in-flight 
transfer alignment and calibration procedure is usually performed prior to release 
of the guided munition munitions. During this period, velocity differences (or 
related quantities) between the aircraft INS and the weapon INS are processed 
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by a Kalman filter to initialize the attitude and heading of the weapon INS, and 
to estimate some of its inertial sensor errors, typically including the gyro and 
accelerometer bias error. The vibration spectrum present at the weapon inertial 
measurement unit is strongly driven by its captive-carry environment, such as a 

5 wing-store station or weapon bay. During the captive carry, the weapon's 
inertial sensor bias errors will be affected by the vibration environment. 
However, after release of the weapon, a significantly different vibration 
environment will exist, leading to gyro and accelerometer bias shifts due to a 
change in the vibration-induced errors. Depending on its magnitude, this change 

10 can negate any benefit of the pre-launch weapon inertial measurement unit 

calibration. There is a need to minimize the detrimental effects of such changes 
in vibration environment. 

Other factors can also cause such bias shifts. Anything that causes a 
change in the vibration environment of a sensor can result in a bias shift. 

1 5 Changes in speed of any vehicle having such sensors can result in different 
vibration environments due to aerodynamic differences and also due to 
mechanically based resonant frequency changes. This can by observed by any 
passenger in a car, especially older cars, as it approaches speeds that are higher 
than intended, or equipment of the vehicle changes through damage. Vibration 

20 environments can also change as engines are operated at different revolutions per 
minute, or depending on the overall load condition of the vehicle. 

There is a need for a better mechanism to compensate for vibration 
induced errors affecting inertial sensors. There is a further need to compensate 
for changes in vibration environments that affect such sensors. There is a need 

25 for such compensation to be performed both accurately and quickly. 

Summary of the Invention 



30 A factory calibration procedure is performed by subjecting a sensor to 

suitable vibration stimulus, thereby exposing sensitivity coeffitifiDlsjelating 
vibration power levels at various frequencies to the biasrlike_error of the sensor. 
In the real-time application of the sensor, aj^^m^^^I^ 
andusfidjfixhaiac^ actual vibrational^ the 
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amount of vibration energy in various portions of the frequency spectrum. 
Predicted sensor v ibration-induced error is calculated based on the measured 
vibration spectrum and the previously- obtained fa ctory calibration^ 
The predicted vibration-induced error is then applied as a compenjation^ 
5 sensor output. 

In one embodiment, high rate vibration measurement samples from the 
accelerometers are provided as the input to a power spectral density estimator. 
The vibration pow er spectrum estimates are then combined with previo usly 
o btained calibration coefficients to geiieTatejrea^ Vibration 
1 0 biassensitiyity coefficjenlsthat are estimated in a contrpjled ^^n^ronmjgntjg 
ensure their accuracy to compensate^ 

In a further embodiment the high data rate vibration samples are ofetaineci 

using acceleration measurement devices such as accelerometers mounted to a 
* ' ' — 

rigid sensor assembly shared with the sensor to be compensated. Vibration bias 
1 5 gstimatesare then combined with outputs from both a plurality (typically three 

or more) of gyros and a plurality .of accelerometers . which combined form an 

inertial measurement unit. 

In yet a further embodiment, a digital signal processor is used to perform 

calculations to obtain the vibration power spectral density and then provide the 
20 sensor vibration-induced error estimates in real time at a high rate. Still further 

embodiments utilize software to perform the calculations, or further 

combinations of hardware, software and firmware. 



Description of the Figures 

25 Figure 1 is a block diagram of a mechanism for estimating and correcting a 

sensor for vibration induced error. 
Figure 2 is a graph showing a vibration-induced error model for a 

simplified example. 

Figure 3 is a flowchart representing a computer implementable method for 

30 estimating and correcting a sensor for vibration induced error. 

Figure 4 is a block diagram showing the estimation of coefficients given 

reference accelerations. 
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Detailed Description 

In the following description, reference is made to the accompanying 
drawings which form a part hereof, and in which is shown by way of illustration 
specific embodiments in which the invention may be practiced. These 

5 embodiments are described in sufficient detail to enable those skilled in the art to 
practice the invention, and it is to be understood that other embodiments may be 
utilized and that structural, logical and electrical changes may be made without 
departing from the scope of the present invention. The following description is, 
therefore, not to be taken in a limited sense, and the scope of the present 

10 invention is defined by the appended claims. 

The description is divided into multiple sections. The first section 
describes a mechanism for estimating and correcting a sens or for vibration, 
i nduced en s*. This is followed by a description of an algorithm implemented by 
the mechanism. A conclusion follows, which describes benefits. 

15 In Figure 1 , a block diagram showing a mechanism forj^timati^g,and^ 

correcrtingifi^ 110. A housing 

or other mounting device 120 in one embodiment contains a micro- 
electromechanical sensor (MEMS) device 124, for example, an inertial 
measurement unit (IMU). This tYPJ^Qi^n^^ 

20 for rate/acceleratio n seeing. The sensor is a Coriolis-type vibratory gyro in one 

embodiment which provides a first output representative of directional changes. ^ \ „ ^ 
An3g<^lerometera26is^ 

repre sentati ve ofaccelgjation. The second output is provided to a vibration 
power spectral density (PSD) estimatorl30, which pjgvidesest^^ 
25 power spectral densities of vibrations as a function of high frequ ency 

ac celeromet er 126 measurementsinSTone embodiment, the frequency of such 
measurements is at least twice the highest expected frequency of vibration to 
ensure sufficient accuracy of measurements. Sample rates may be on the order 
of 5KHz or higher. 

30 The PSD estimates are formed at a relatively low rate, for example 1 Hz, 

and are provided to a compensaUon ^ 

coefficients IJiLto. provides .>3bration^duced bi^ estimates or compensation 
values , both for the gyroscope first output and the ^££|£^^^f^S^^I^!Uv 
An adder 150 receives the first output, and a gyro vibration-induced bias 
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estimate and provides a compensated gyroscope output on a line 155. An adder 
160 receives the second ouipu i_and a fl,.fl£££ljgro^^ * 
estimate from modu le J 3,4^md outpuUm 
a jine 16.5^ 

5 To ensure that the vibration measurements correspond to the sensor being 

compensated, the components are mounted to a rigid sensor assembly indicated 
at 1 70 and shared with the sensor to be compensated. 

The basic model for vibration-induced bias error is that the bias depends 
in part on the vibration power spectral densities (and possibly cross power 
10 densities): 

vibration bias = f(Sxx (co), Syy (co), Szz(co)) 
where S xx (co), Syy (co), and Szz(co) represent the vibration power spectral 
densities in the sensor's x, y and z axes. In one embodiment, these accelerometer 
axes are those that are customarily referred to as the input axis, output axis, and 

1 5 pendulous axis of the accelerometer. 

In one embodiment, an accelerometer exhibits a sensitivity to input axis 
vibration, a common characteristic of devices with non-negligible scale factor 
nonlinearity. To simplify, assume that the device does not have cross-axis 
vibration sensitivity. In this case, the vibration induced bias 6a x depends only on 

20 input axis vibration: 

5a x =f(S xx (<3)) 

The vibration spectrum can be broken up into a set of frequency bins, and the 
sensitivity of the bias Error to the vibration power within each bin is defined as 
illustrated in Figure 2, which is intended as a typical graph of the spectral density 
25 versus frequency. The vibration-induced bias is then expressed as a sum over 
the frequency bins: 

6a x =Zi Ci& 2 

where gj 2 represents the vibration power (mean squared acceleration) in the i tb 
frequency bin, and c< is the vibration-induced bias error sensitivity coefficient for 
30 that bin. Estimates for the coefficients Cj shown in the model are determined 

* ■ 1 11 ■t" ■' ■■■■■ ■i f • ■r—rTwr MwyW rwa 

d uring &g ^2iyf^^!^- Q " °^ SlSjS^S^dJBS^jgj sensors. 

In an alternative embodiment, the coefficients are determin ed for a small 
number of sensors referred to as a family of sensors, and the typical coefficient 
values are applied as compensat ion coeffi cients for MLs ensors of th e same type. 
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The calibration process may thus be applied to a single manufacturing batch of 
sensors, or any other family quantum of sensors that are likely to be similarly 
affected by vibration. This eliminates the need to perform extensive factory 
calibration procedures on every production sensor, but may retain most of the 

5 accuracy of the vibration-induced error compensation. 

In one embodiment, the frequency bins are of equal size and cover the 
frequency range at which vibration power may significantly affect the sensors. 
In further embodiments, the frequency span of the bin is adjusted such that in 
frequency ranges where vibrations are most likely to affect the sensors, more 

1 0 bins are concentrated, and in other less significant frequency ranges, fewer bins 
are used. In one embodiment, the bins are coupled to a bank of bandpass filters, 
each covering one the frequency bins. The mean squared acceleration at the 
output of each of the bandpass filters is used as a measure of the vibration power 
in each bin. 

15 a sine sweep vibration profile is used as test stimulus to provide 

observability of the model coefficient s. The center frequency of the narrow- 
band sinusoidal vibration stimulus is swept over the range of vibration frequency 
typically encountered by the sensor in its intended application. 

A block diagram flow chart is shown in Figure 3. Measured 
20 accelerations (dv/dt) for x, y and z axes are provided at 310 and sampled in each 
of the respective axes at 312, 314 and 316 over a desired time interval. 
Respective power spectral densities are then computed at 322, 324 and 326. For 
each axis, the vibration power is then computed for each desired bin at 332, 334 
and 336. There are N bins shown, and the actual number of bins varies 
25 depending upon the desired accuracyof correction and computational power 

available. The vibration induced errors are calculated at 340 with predetermined 
vibration sensitivity coefrid ents^provided at 345. The calculation of the 
vibration induced errors at 340 is discussed in further detail below, and provides 
the coefficients for each of the axes at 350. 
30 Figure 4 is a block diagram showing calibrating of a vibr^^induced 

error model for a tria^fcaceel^ including the dependence of each 

sensor's error as a function of vibration power in a set of three principal axes that 
are nominally aligned with the accelerometer's input axes. Measured 
accelerations in each axis are provided at 410 in response to actual reference or 

6 
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truth accelerations provided by a current test condition at 412. A difference for 
each axis is provided by adders 416, 418 and 420, which are provided to a 
vibration sensitivity coefficient estimator 430. The measured accelerations 410 
are also provided to a buffer 435, which samples the accelerations in each axis 

5 over a predetermined calibration time at a rate at least twice the highest 
frequency desired to be sampled in order to obtain valid sampling data. 
Vibration power spectral densities are then computed at 440 along each axis and 
then the vibration power in each of N frequency bins is computed at 445. 
These values are then provided to the estimator 430. 

1 0 Calibration of the necessary model requires a number of different test 

conditions, where a test condition is defined by the axis of the applied vibration 
stimulus and the frequency bin in which most of the vibration energy is 
concentrated. In this case, a model for the vibration-induced error on the jth test 
condition is 
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which can be written as 
8aj = Hjc 

20 where c is the 9N-element column vector formed form the vibration s ensitivity 
coefficients 

Cxxi, Cxyj, Cxzi, Cyx,, Cyyi, Cyzj, Czx i5 Czy i5 Czzi for i = 1, 2 . . . N 
and Hj is a 3 x 9 N coefficient matrix comprised of the measured 3-axis vibration 
power levels in the N vibration frequency bins. 
25 Data collected over a sufficient number (M) of test conditions (vibration 

stimulus in each of the three principal axes, each over a number (3 N) of 
different dominant vibration frequencies) can be used to develop the composite 
vibration-induced error model 



WO 01/11318 



PCT/US00/21396 









An 
OCL 2 













or more compactly, 
y = Hc 

Here the measurement vector y represents the vibration-induced 
acceleration errors present in the three accelerometer's outputs at each test 
5 condition. The vector of unknowns (c) represents the vibration sensitivity 
coefficients to be estimated. The matrix H is comprised of the measured 
vibration power levels in the three principal axes and the 3 N vibration 
frequency bins. 

Given the form of the above measurement model, the well-known least 
10 squares estimate of the vibration sensitivity coefficients is given by 
□ = (H T Hy l H T y 

Other well-known estimation algorithms are also used to generate the 
vibration sensitivity coeffi cient estimate s in different embodiments, including a 
weighted least squares estimator, maximum likelihood type estimators, and 

1 5 recursive estimators such as the Kalman filter. 

As an alternative to factory calibration, some applications use on-line 
estimation of the vibration sensitivity coefficients. In one embodiment, a global 
position system (GPS) aided inertial navigation system, the GPS provides 
accurate measurements of position and velocity (or equivalently, pseudorange 

20 and/or deltarange for a plurality of GPS satellites). The GPS measurements 

expose the position and velocity errors of the navigation solution formed using 
the IMU gyroscopes and accelerometers, which are affected by host vehicle 
vibration. Processing the GPS measurements and the measured vibration 
spectrum in, for example, a GPS/inertial Kalman filter algorithm provides 

25 observability of the vibration-induced error sensitivity coefficients. 

Although the simplified example assumes that the sensor is affected by 
only a single component of the three dimensional vibration stimulus, in practice 

8 
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it is expected that a sensor may be influenced by the vibration spectra acting on 
all three axes. This three-dimensional vibration stimulus would be characterized 
by the vibration PSDs in the sensor's three principal coordinate axes, and by the 
associated cross power spectra. 

5 

Conclusion 

It is to be understood that the above description is intended to be 
illustrative, and not restrictive. Many other embodiments will be apparent to 
those of skill in the art upon reviewing the above description. For instance, 

10 while inertial sensor compensation is described, the outputs of virtually any type 
of sensor, such as pressure, gas flow, temperature, etc. may also be compensated 
in the same manner, provided that suitable measurements are available to 
characterize the real time vibration environment. 

In one embodiment of the invention the vibration power spectra 

1 5 estimation function is performed using software algorithms hosted in a general 
purpose microprocessor. An alternative embodiment uses a dedicated digital 
signal processor type hardware device as a coprocessor tasked with generating 
the PSD estimates. Still further embodiments will be apparent. 

In further embodiments, sets of calibration coefficients are provided for 

20 selected temperature ranges for sensors which exhibit temperature-dependent 

vibration sensitivity coefficients. This embodiment includes a temperature sensor 
located proximate to the sensor being compensated. An output of the 
temperature sensor is used by the compensation module to obtain the correct set 
of compensation coefficients. 

25 The scope of the invention should, therefore, be determined with 

reference to the appended claims, along with the fiill scope of equivalents to 
which such claims are entitled. 
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What is claimed is: 

1. A method of reducing the effect of vibration-induced errors in sensors, 
5 the method comprising: 

obtaining high data rate vibration samples (126); 
estimating vibration power spectral density (130); and 
combining the estimates with previously determined calibration 
coefficients to generate real time vibration-induced error correction (134). 

10 

2. The method of claim 1 wherein estimating vibration power spectral 
density further comprises computing vibration power values in multiple 
frequency bins. 

15 3. The method of claim 1 wherein the predetermined vibration sensitivity 
coefficients are determined in a factory vibration stimulus calibration procedure. 

4. The method of claim 1 wherein the previously determined calibration 
coefficients are determined with respect to the actual sensors being corrected. 

20 

5. The method of claim 1 wherein the previously determined calibration 
coefficients are determined with respect to a family of sensors. 

6. The method of claim 5 wherein the family of sensors comprises sensors 
25 manufactured in a same batch. 

7. A computer readable medium having instructions stored thereon for 
causing a computer to perform the method of claim 1. 

30 8. The method of claim 1 wherein the power spectral density is measured in 
one axis. 

9. The method of claim 1 wherein the power spectral density is measured in 
three axes. 
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10. A computer implemented method of compensating for vibration induced 
error in sensors, the method comprising: 

buffering vibration samples in at least one axis of vibration; 
5 computing a power spectral density for such buffered vibration samples; 

computing vibration power values in multiple frequency bins; 

combining the vibration power values with predetermined vibration 
sensitivity coefficients to obtain estimated vibration induced errors. 

10 11. The method of claim 1 0 wherein the frequency bins are equal in 
frequency range. 

12. The method of claim 10 wherein the frequency bins are smaller in 
frequency range about a predetermined frequency. 

15 

13. The method of claim 10 applied to multiple axes of vibration. 

14. The method of claim 1 0 wherein the power values in the multiple 
frequency bins are representative of the mean-squared vibration level in each 

20 bin. 



1 5. The method of claim 10 wherein the vibration samples are obtained via 
band pass filters corresponding to the frequency bins. 

25 1 6. The method of claim 1 0 wherein the predetermined vibration sensitivity 
coefficients are determined in a factory vibration stimulus calibration procedure. 

17. The method of claim 10 wherein the predetermined vibration sensitivity 
coefficients are determined via comparison of measurements unaffected by 
30 vibration to those obtained from the sensor which is affected by vibration in real 
time. 
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18. The method of claim 1 7 wherein the sensors comprise a gyroscope and 

accelerometer, and wherein the measurements unaffected by vibration comprise 
T — - 1 ,, W M m > 

GPS position measurements. 

19. A computer readable medium having instructions stored thereon for 
causing a computer system to perform the method of claim 10. 

20. A method of calculating correction coefficients for an inertial sensor 
comprising: 

measuring acceleration of the sensor; 

computing vibration power in selected frequency bins corresponding to 

the measured acceleration; 

providing a difference in measured acceleration from a reference; and 
combining the measured acceleration and the difference to obtain 

vibration sensitivity coefficients. 

21 . A system that provides vibration calibration coefficients for a sensor 
affected by vibration, the system comprising: 

means for receiving measured acceleration; 
means for computing vibration power in selected frequency bins; 
means for providing a difference in measured acceleration from a 
reference; and 

means for combining the measured acceleration and the difference to 
obtain vibration sensitivity coefficients. 

22. The system of claim 21 and further comprising means for ensuring that 
the sensor is subjected to the measured acceleration. 

23. A system for compensating for vibration induced sensor (120) error, the 
system comprising: 

an accelerometer (126); 

a power spectral density estimator (130) coupled to the accelerometer 

(126); 
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a compensation module (134) coupled to the power spectral density 
estimator ( 1 30) and to a plurality of calibration coefficients (140) that provides a 
compensation value for a sensor output. 

5 23. The system of claim 22 wherein the accelerometer comprises an 
accelerometer for each axis of vibration. 

24. The system of claim 22 wherein the sensor comprises a vibratory 
mechanism. 

10 

25. The system of claim 24 wherein the sensor comprises a gyros. 

26. The system of claim 25 wherein the sensor further comprises an 
accelerometer for which compensation values are also provided by the 

1 5 compensation module. 

27. The system of claim 22 wherein the sensor comprises an inertial sensor. 

28. The system of claim 27 and further comprising a rigid assembly for 
20 mounting the inertial sensor and the accelerometer. 

29. A system for compensating for vibration induced sensor error, the system 
comprising: 

means for buffering vibration samples in at least one axis of vibration; 
25 means for computing a power spectral density for such buffered vibration 

samples; 

means for computing vibration power values in multiple frequency bins; 
means for combining the vibration power values with predetermined vibration 
sensitivity coefficients to obtain estimated vibration induced errors. 
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